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The  Average  First  Recurrence  Time 


Introduction 

The  kinetic  theory  of  gases  contains  a  fundamental 
paradox  in  that  the  results  obtained  from  the  theory  are  ir¬ 
reversible  despite  the  fact  that  gases  are  assumed  to  be  com 


posed  of  molecules  which  obey  the  reversible  laws  of  mechanics. 
The  paradox  is  further  heightened  by  Poincare’s  Cycle  Theorem 
which  states  roughly  that,  no  matter  what  is  the  initial  con¬ 
figuration  of  the  molecules,  the  gas  must  return  to  a  state 
which  is  arbitrarily  close  to  the  initial  state  and  that  it  must 
do  so  infinitely  often.  Physicists  have  tried  to  explain  this 
paradox  by  stating  that  the  time  for  one  such  Poincare  cycle  is 
so  large  that  the  recurrence  of  an  Initially  improbable  state  is 
unlikely  during  the  times  normally  available  for  observation. 

For  a  fuller  discussion  of  these  topics  we  refer  to  an  article 
by  Chandrasekhar1. 

In  the  present  paper  we  consider  the  first  recurrence 
time  of  a  very  3imple  dynamical  system,  namely,  one  which  has 


k+1  degrees  of  freedom  x-^Xg* 
periodic  in  the  time  t  .  Put 


,xk+l' 


each  of  which  is  simply- 


i(  T  ) 


^2iriv  .t 

J  > 


3=1,2, ... ,  k+1  j 


then  the  problem  is  to  determine  the  smallest  positive  value  of 
t  such  that  the  point  (x^(x  )j*»*jXk+1(x  ))  will  be  within  a 
preassigned  neighborhood  of  the  initial  point  (x1(0) , . . .  ,xlc+1(0) ) . 
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For  simplicity,  the  neighborhood  is  assumed  to  be  such  that 


(1) 


lvjT"mjl  <  e'  ^1*2,..., k+1, 


where  are  integers.  The  smallest  positive  value  of 

r  such  that  the  inequalities  (1)  are  satisfied  for  some  set  of 
integers  will  be  called  the  first  recurrence  time  of  the  system, 
Note  that 

T  =  T  (Vt . Vwi  £)• 


We  wish  to  find  the  average  of  the  first  recurrence  time  for  all 
such  systems,  that  is, 

1  1  , 

(2)  Tav  «  f  ...  f  t( v]/  •  *  •  *  vk+l^  €  )^vi •  •  ,(^vic+l 

o  0 

and  in  particular  the  asymptotic  behavior  of  t  „  as  e  ap- 

a  V 

proaches  zero. 

Suppose  that  is  the  largest  of  the  and  instead 

of  the  inequalities  (1)  consider  the  following: 

(3)  .  Vk+1T ”mk+l  "  °'  iVjT~mjl<e'  J=l>2,...,k. 

Put  <jj  «  vj/vk+l'  mk+l  “  then  (3)  becomes 

(4)  |  tOj-m^  |  <  e,  j*=l,2, . . . ,k, 

where  t,  are  integers.  We  put 

/ 

i  1  . 

tav  *  /  /  t(a1, . . .,ak;  ejdc^ . . *dak. 

o  o 
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We  shall  show  that  there  exist  two  constants  c1  and  c2  such 
that 


(5) 


cle"k  <  ‘av  <  °2e'k 


For  the  case  k*l,  more  can  be  obtained.  We  prove  then 


that 


(6) 


tgv  -  log  2<=  _1  +  0(€  -1/2). 


The  Qne-Dlitienslonal  Case 

For  any  a  in  the  unit  interval,  define  t  «  t(a>  6  )  as 
the  smallest  positive  integer  such  that  |  ta-m  j  <  e  for  some 
integer  m,  i.e.,  such  that  ta  is  within  €  of  an  integer.  Our 
problem  is  to  evaluate 


(7) 


/  t(a, €  )da. 


Let  the  integer  n  be  chosen  so  that  —  <  e  <  -1—  .  Then  it  is 

n  “  n+1 

clear  that  for  any  a. 


t(a,  i)  <  t(a,  €  )  <  t (a, 


n  = 


and  so 


1  l  1  1  l 

/  t(a,  i)da  <  /  t(a,  e)da£/  t(a,  ~-)dd. 


n’ 


Thus  we  examine 


/  t(a,  i)da. 


I  * 


-  4  - 


Any  real  number  a  is  said  to  have  an  admissible  value  k 
(a  positive  integer)  if  ka  is  within  1/n  of  an  integer.  Thus 
t(a,  1/n)  is  the  minimum  of  the  admissible  values  for  a*  Any 
rational  number  a/b  with  0  <  a/b  <  1  defines  an  interval 


l(aA>) 


®+  i.) 

b  nb' 


y 


and  all  points  in  l(a/b)  have  admissible  value  b.  Furthermore 
b  is  an  admissible  value  for  no  real  values  of  a  other  than 
those  in  the  intervals  l(0/b),  l(l/b),  l(2/b) , . . . ,l(b/b) . 

Now  if  t  is  the  minimum  integer  such  that  ta  is  within 
1/n  of  an  integer  m,  then  ta-m  =  ^  with  j  5  |  <  1.  Thus 

a-m/t  «  J-  ,  and  m/t  is  in  its  lowest  terms  since  otherwise 
tn 

t/(t,m)  would  be  an  admissible  value  for  a  smaller  than  t  itself. 
Thus  a  e  I (m/t). 

Let  Fn  denote  the  Farey  series  of  order  n. 


*.01  n-1  1 

n  1  n  n  1 

the  series  of  rational  numbers  a/b  in  ascending  order  satisfying 
0  <  a  <  b  <  n  with  (a,b)  *»  1. 

Now  if  the  Farey  series  of  fractions  a/b  say,  is 

used  to  define  a  collection  C  of  intervals  l(a/b),  we  shall 
prove  that  C  covers  the  unit  interval,  but  that  no  single 
interval  I  (a/b)  contains  any  member  of  other  than  a/b.  From 

this  it  follows  that  to  find  t(a,  1/n)  we  need  merely  locate  a 


al  a2 

between  two  consecutive  fractions  ?-=■  and  —  of  Pw  ,  and  we 

b2  n-1 

have 

3  3 

t (a,  l/n)  =  b,  if  a  is  in  I(~)  but  not  in  I(~); 

bl  2 

ap  a, 

t (a,  l/n)  =  bp  if  a  is  in  I(r~)  but  not  in  I(_±); 

b2  bl 

3  3 

t(a,  l/n)  =  min(b1,b2)  if  a  is  in  both  I(~)  and  l(-£). 


To  prove  the  assertions  above  about  the  collection  of 
al  a2 

intervals  C:  If  —  and  —  are  adjacent  fractions  in  F„  , 

Vi_  b.  n-1 


al  a2 

then  ~  ~  «  + 

h  bz  " 


h  b  n-i 

bl  d2 

1  a2 
-  — -  and  since  b~  <  n  then  r-—  is  not  in 

'  b1b2  d  b2 


3  3  3 

l(rr:-).  To  prove  that  I(t~)  and  I(r^-)  meet,  (which  will 
D1  D1 

establish  that  C  covers  the  unit  interval)  we  must  prove  that 
— +  -i—  >  ,  This  follows  from  the  fact  that 

nb1  nbg  =*  b-jbg 


bi  +  bg  >  n. 


al  a2, 


From  this  argument  we  see  that  the  interval  (— *-— ),  of 

bl  D2 

length  — — ,  splits  into  two  parts  as  regards  the  value  of 
blb2 

t(a,  i): 

t(a,  i)  »  min(b, ,bp)  over  a  subinterval  of  length  - - - 

n  n.min(b1,b2) 

t(a,  — )  »  max(b,,b0)  over  one  of  length  — - - - - -  . 

h  12  b^bg  n*min(b1,bg) 


min(b1Jb2) 

n*min(b1Jb2) 


+ 


maxCb-^bg) 

n.min(b1,b2) 


a2/^b2 

/  t(ct,  i)da  » 
n 


maxCb-^bg) 


blb2 


n-max(b^,bg) 

n*min(b1Jb2) 


Let  us  denote  minfb-^bg)  by  b  and  maxCb-^bg)  by  B,  so  that 


(8) 


/  t(a,  ^)da 
o  n 


n-B\ 
nb  ' 


where  the  sum  ranges  over  all  pairs  of  consecutive  fractions  in 
Fn-1‘ 

Lemma  1.  Consider  the  set  SK  of  all  fractions  in  F  , 

d  n-i 

with  fixed  denominator  b,  where  1  <  b  <  n-2.  Let  T  be  the  set  of 

ss  c? 

all  fractions  which  are  neighbors  to  members  of  Sb,  and  which 
have  denominators  exceeding  b.  If  B  is  any  integer  such  that 
(B,b)  «=  1 ,  b  +  B  >  n,  n-1  £  B  >  b  +  1,  then  the  set  T  has 
exactly  two  members  with  denominator  B. 

Proof.  The  equation  xb  -  yB  -  1  has  exactly  one  solution 
in  integers  x  and  y  such,  that  1  ^  x  <  B-l,  1  £  y  <  b-1.  The 
fractions  x/B  and  y/b  are  in  ?n-1>  and  are  neighbors,  for  if 
y/b  had  the  neighbor  h/k  on  the  right,  y/b  <  h/k  <  x/B  then 
hb  -  yk  «  1  and  xb  -  yB  **  1  imply  that 

b(xk  -  hB)  «*  k  -  B. 


dL; 


1 


t 
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But  xk  -  hB  >  1.  so  that  k  -  B  >  b  and  k  >  b  +  B  >  n,  a  contra- 

S3  22  S3  SC 

diction. 

Similarly  the  equation  Xb  -  YB  <=  -1  haa  exactly  one 
solution  X,  Y  with  1  <  X  <  B-l,  1  <  Y  <  b-1  and  likewise  the 
fractions  X/B  <  Y/b  are  neighbors  in  F  Thus  we  have  shown 

that  the  set  T  contains  at  least  two  fractions  with  denominator 
B.  On  the  other  hand  if  a/b  and  A/B  are  neighbors  in  Fn 
then  aB  -  Ab  *»  +  1,  so  that  the  two  cases  we  have  examined  are 
the  only  possible  ones.  Hence  the  lemma  is  proved. 

Returning  to  the  formula  (8),  we  first  compute  the  sum 
S  1/n.  Let  v'  (k)  be  the  Euler  function  denoting  the  number  of 
integers  less  than  and  prime  to  k,  and  define 


<|(n)  "  <f(l)  +  f  (2)+*  •  •+  f  (n)  • 

Then  the  number  of  terms  in  Pn-1  with  denominator  k  <  n-1  is 
<j 9  (k),  except  in  case  k  «  1,  there  being  two  terms  with 
denominator  1.  Hence  we  have 


S  n  ”  n  [  +  <^2)+,,,+  ffa"1)  j  -  £  I  (n-1). 

By  virtue  of  this  and  Lemma  1  we  see  that  (8)  can  be  written  as 


/  t(a,  i)da 

% 


if  (n-1)  +  i 
n  v  '  n 


(9) 


n-2 

2 

b=»l 


H'1  n-B 

B**b+1  b 
B+b  >  n 

(B,b)  -  1 


H 


i s 
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•zm^ 

Lemma  2.  (F(m)  33  — r—  +  0(m  log  m) 

x  TT^ 

+  J£l£l  +. . .+  ^~tm)  »  +  O(log2m) 

12  m  71-2 

JfJg+11  +  (m+2)  +  . . .+  (,2m).  „  £j°S_2  +  0(logJR)  . 

(m+1)2  (m+2)2  (2m)2  ir2  m 

2 

Proof.  The  first  Is  in  Hardy  and  W.  ight  .  The  third 

3 

result  Is  a  special  case  of  a  result  of  D.  H.  Lehmer  .  The 
second  result  can  be  derived  from  the  first  as  follows: 


+  X121.+...+ 

1  a  m 

_  1(1)  |  1(2)-  1(1)  |  1(3)- £(a)  +...+  Km)-I(m-l) 

12  3  m 

-"z1  g(f),  +  JlBl 

>1  J(J+1)  m 

-  % ("'i1  ♦  2!  (+  o(5LiHS-B)  +  ^  0(ii2iLi) 

TT2  [  J-l  j(j+l)  m  J  m  J=1  J(J+1) 

■  ~  /  2m  +  O(log  m)^1  +  O(log  m)  +  0(log2  m) 


%  +  0(log2  m). 
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Lemma  3.  Let  (f  ( j,n)  denote  the  number  of  integers  <  J 
that  are  prime  to  n.  Then 

(?(3>n)  *»  J  +  0( V”1n)  . 

Proof.  Lehmer  ![3,  p.  1188}  proves  that 
|  ((  (J>n)  -  j  |  <  d(n) 

where  d(n)  is  the  number  of  divisors  of  n.  To  prove  that 
d(n)  »  0(\Tn),  say  d(n)  <  4v/n  to  be  definite,  we  can  use 
induction  on  the  number  of  prime  factors  of  n.  Write 

kl  ks  2  kl 

n  *»  p^  •  •  •  pg  •  If  s  a  1,  we  must  prove  that  (l+k.^)  <  ^p^^ 

It  can  be  readily  verified  that 
(*)  (H-l^)2  <  4Pl  1 

for  p  =«  2,  and  so  a  fortiori  for  p  >  2.  If  s  a  2,  we  must  prove 
that 

(l+k1)2(l+k2)2  < 

This  is  a  consequence  of  (^  ).  Next,  if  we  assume  the  result 
for  numbers  with  fewer  than  s  prime  factors,  with  s  >  3,  then 

s 

we  have 

d(n)  -  (l+ItgJdf-Sj-)  <  (l+k3)  4  / — — —  <  4/n 

Pa  3  ‘Pa  3 

2  ks 

provided  that  (l+ks)  <  pg  ,  and  this  is  readily  verified 
since  p_  >  5. 

rS  » 


s  10 


r 


Lemma  4.  Let  S(j,n)  denote  the  sum  of  the  integers  that 
are  <  J  and  prime  to  n.  Then 

2 

S(j,n)  -  J— f'M  +  0{j  fn). 

2n 

Proof.  We  use  induction  on  the  number  of  prime  factors 
of  n.  Let  us  write 

k,  k 

n  «  Pi  •  •  *P8  ,  n!  -  — Y"  '  -  QPS  +  0  <  r  <  p3 . 

Ps  3 

First,  we  prove  the  identity 

(10)  S(j,n)  ~  S(j,n1)  -  PgSCq,^). 

This  can  be  seen  by  observing  that  none  of  the  numbers 
p  _ ,  2p  . ...,qp_  are  included  in  S(j,n),  but  some  of  these  are 
included  in  SO,^),  namely,  the  Integers  in  l,2,***,q  that  are 
prime  to  n^,  each  multiplied  by  pg . 

If  s  ■  1,  then  (10)  gives  us 

S(J.n)  -  S(J,1)  -  p3s(q,l)  -  -Pial|±il  . 

Now,  since 

2n  2  Pi 

2 

S(j,n)  -  »  ~(J  -  Pi^'PiS  +  J2/Pi) 

2n  2 

-rq  +  l+  l  *^-<pq  +  |  +  |*pq  +  r«J  <  J  fn  . 


we  have 


Before  turning  to  the  rest  of  the  induction  argument,  we 

handle  the  special  case  of  Integers  n  divisible  by  2  and  3  but 

by  no  other  primes .  Thus  *»  •  A  ,  and 

n  o  3 


S(J,n)  »  S(j,6)  **1  +  5+  7  +  11  +  13+  17+  •  •  • 

«  6  +  18  +  3CM-***+[  2jlj  -lj  6  +  0(j) 

-  6^1  +  3  +  5+-*-+(2 /s]  "1)J  +  °(j) 

-  6/i f  +  0(J)  -£  +  0(J)  -  +0(jVTi). 

lO-i  6  9n 


We  now  use  induction,  assuming  that  the  lemma  holds  for 
integers  n  with  fewer  than  s  prime  factors;  and  we  may  presume 
that  ps  >  5  because  of  the  above  special  case.  Applying  (10)  we 
have 


S(J,n)  «  SCJ,^)  -  PgSCq,^) 

-  -  P°q  y’(ni)  +  0(  j  fnx)  +  ps0(q 


2n, 


2n, 


But 


lllnl  .  illUll  (1_  I_)  .  itlbl  -  (,2p.+  2qr+i 


2n  2n, 


Ps  2n, 


Hence 


S(J,n)  -  illSl  .  Ipl  (qr  +  JL)  +  OU/T.)  +  0(p  q  /n^  ) . 
2n  nl  2p  A  8 


_  *  ji-) 
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(ii) 


Ti  - 


n-2 

b-k 


n-l 

2 

B=b+1 

(Bjb)=l 


n-B 

b 


In  T2  we  have  2b  <  n,  n  >  2b-l,  so  B  +  b  >  n  implies 
B>n-b>b+l;  therefore  the  sum  ranges  from  B  =  n-b  to 

B  =  n-l  and  we  write 

f  h-1] 

LtJ  n-l 


(12)  T2  =  2 


n-B 

b=l  B=n-b  b 

(B,b)r=l 


We  begin  with  (11).  Replace  B  by  b  +  k  in  the  inner  sum;  then 


n-l 

2  n— B 

B=b+1  b 

(B,b)-1 


n 

2  n-B 

B=b+1  h 

(B,b)-1 


n-b-k 


n-b 

2  . 

k=l  b 

(k,b)=l 


y  n~b  v  i£ 
b  b 


2=^  <j>  (n-b,b)  -  i  S(n-b,b) 

S±  f  ISdilJM  +  o(fb)j 

<p(b)  +  0(f n) 


-  1  fin-bfjM  +  o[(n-b)fb||+o(l) 


2be 


with  the  help  of  Lemmas  3  and  4. 


Therefore,  using  Lemma  2,  we  find  that 


T  "j2  (£  <£M.n£M+  <n*l+0(rn)l 
1  b.[2±i]1'2  b2  b  2  3 

L12L1  -  n(%  -  3|)  ♦  I(2|i  -  52!)  +  O(nrn) 

2  IT2  IT2  7T2  2  7T2  4tt2 


^  (3  log  2  -  ^)  +  0(n  fn) 


14 


To  evaluate  (12),  we  begin  with  the  inner  sum.  Define  r 
as  the  least  positive  residue  of  n  modulo  b,  thus 
r  £  n(mod  b),  1  <  r  <  b.  Then  we  replace  B  by  n-r-i  so  that 
(B,b)  =  1  becomes  (i,b)  =  1. 


We  have 

H"1  n-B 

B=n*-b  b 

(B,b)=l 


r+i 


b-r 
=  2 

i«  -r+1  b 
(i,b)=l 


2  |  +  2  ~ 
b  b 


(13) 


b-r 

b  i-l 

(i^b)-l 


|  +  1 

b  i=»  -r+1 


<1  +  |) 


-  (f  (r-l,  b). 


(l,b)-l 


But,  since 


b-r 


-1  ,  b-r  ,  b-1  i 

2  i  +  2  (1  +  .i)  -  Z  i  +  2  V 

i-l  b  i=  -r+1  b  i-l  b  i=b-r+l 

(i,b)-l  (i,b)-l  (i,b)=l  (i,b)«l 


b-i 

!  z  £ 
1  b 

(i,b)-l 


Vf\  .  tfV, 


S  % 


+  2  . 

1  u  i=l  b 

(i,b)-l  (i,b)-i 


.  1M. , 

2 

we  find  that  (13)  becomes 

r  O’).  +  (r,b)  +  0(1) 

b  2 

-  im  +  0(rb) 


by  Lemma  3. 


2 


We  use  this  result  to  evaluate  the  sum  Tg*  We  have 
n-1 


m 


T2  ”  A  I 


+  0(  rb) 

i  ^ 


b-1  l  2 
n>2 


2  2 
*  tr 


+  0(n^  n)  by  Lemma  2 


Then 


3n  \ 
^2  +  °(n  /n) . 


si  L 


3  ) 


T1  +  t2  -  A  f3  106  2  -  fj  +  0(n/n) 
Thus  the  equation  (9)  becomes 


/*<“'  i)da  -if  (n-X)  +  £  [3  log  2  -  |  £  +  o(  /n) 


( 


n’j  ^T+  °(n  106  +  7a{6  log  2-s^+  0(T5) 

+0(1-5). 


7T 


Consequently,  we  have 


/  t(o,  e  )da  -  — -f  2  (~)  +  0(  c  “1//2) 

)  ir2  * 


which  proves  ( 6  ) • 
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To  illustrate  the  higher  dimensional  cases,  we  take  the 
two  dimensional  situation.  The  integer  b  is  an  admissible  value 
for  all  pairs ' of  real  numbers  lying  in  the  region  defined  by  the 
Cartesian  product  of  any  interval 


l(0/b) ,  l(l/b),...,l(b/b) 


with  any  other.  These  intervals  have  length  2/n  (disregarding 
what  falls  outside  the  unit  interval  of  course)  and  the  region 

p 

defined  by  the  Cartesian  product  has  area  4/n  .  Thus 


1  1 

/  /  t  (a-,  ,cXp  t 

o  o 


i)d0lda2  >,1^+2 


a/  4 


8 


because  the  region  with  actual  (not  admissible)  value  b  is 

4 

at  most  of  area  .  Thus  the  t-values  larger  than  -f-  have 
-  n2  4 

been  replaced  by  smaller  values.  In  k  dimensions,  this  argument 
proves  the  existence  of  a  positive  constant  c,  such  that 

1  1  1  1  t  k 

/  /  /  •**  /  t(a1,...,ak,  — )da^  •••  dak  >  cxn  . 

o  o  o  o  n 

Ir 

This  integral  is  less  than  n  ,  because  by  the  pigeon  hole 

If 

procedure  we  can  find  a  positive  integer  m  £  n  for  any 
a1,...,ok  such  that  ma1,...,mak  are  all^  within  —  of  an  integer. 


IT  - 


therefore,  t  <  n  arid  then 


11  1 
/  /  • • •  /  t  da  *  * ‘  aa^  < 


,  _k 

n  , 


O  O 


This  completes  the  proof  of  (5). 
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